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On . Many theories for the formation of massive black holes (BHs) in galactic nuclei 

' predict that the stellar motions should be anisotropic within the radius of in- 

fluence of the BH, r/i = GMh/cr^- If the BH forms slowly by the accumulation 
O , of gas, the stellar action variables are adiabatically conserved and the orbits 

I become gradually more circular (Young 1980; Goodman & Binney 1984). The 

l/^ ' induced anisotropy is quite modest, however, with /? = 1 — cr^ /a^ typically lying 

between —0.3 and in models where the stellar motions are initially isotropic 
(Quinlan et al. 1995). In another class of models, nuclear BHs grow through 
^ ■ the accretion of other BHs acquired from galaxy mergers (e.g. Ebisuzaki et 

'nI" , al. 1991). Dynamical friction drags the BHs into the nucleus where they form a 

bound pair. The binary separation continues to shrink, at first due to three-body 
scattering processes in which the BH binary ejects stars, and later through grav- 
itational radiation (Begelman et al. 1980). Because stars on elongated orbits 
are more likely to be ejected, three-body scattering introduces a circular bias 
in the velocity distribution of the non-ejected stars (Quinlan 1996a). Detailed 
simulations of this process (Quinlan 1996b) suggest that the induced anisotropy 
could be much greater than in the adiabatic growth model. 
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^ • Detecting velocity anisotropy in hot stellar systems is generally a difficult 

c/3 . task, requiring accurate measurements of the line-of-sight velocity distribution 

^ ' as a function of position. However if the functional form of the gravitational 

potential is known, the anisotropy can be inferred from the stellar velocity dis- 
^ i persion profile alone (Binney & Mamon 1982). In the vicinity of the BH, the 

' potential of a spherical galaxy can be approximated as 

= — — + ( xj ^ ^ 

where M/L is the mass-to-light ratio of the stars and <I>l is the "potential" 
corresponding to the stellar luminosity distribution. Eq. (1) is not completely 
general, since it assumes that M/L is independent of radius. However it is likely 
to be accurate within the region where the gravitational force is dominated by 
the BH. M/L can be determined from the velocity dispersion profile via the 
virial theorem, independent of any assumptions about the velocity anisotropy. 

Based on observations of a 1" ionized gas disk. Ford et al. (1994) and Harms 
et al. (1994) inferred the presence of a 2.4±0.7 x 10^ Mq BH in M87. The radius 
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of influence of this BH would be 80 pc f« 1"; but it would dominate the 
gravitational force out to a distance of roughly Sr^. Thus we expect Eq. (1) to 
be accurate within the region where the BH could affect the stellar kinematics. 

Fig. 1 shows estimates of l3{r) for the stars in M87 based on van der Marel's 
(1994) velocity dispersion measurements. The data were corrected for seeing and 
for instrumental blurring using a regularized algorithm (Merritt & Oh 1996). For 
Mh > 1.0 X 1O^M0, the stellar motions are significantly anisotropic, at > ar, 
within 1" — 2". The anisotropy is greater than predicted by the adiabatic growth 
models but may be consistent with the predictions of the binary BH model. 
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Figure 1. Anisotropy parameter (3 = 1 — /crj; as a function of radius for tlie stellar motions 
in M87. Curves are for assumed black hole masses M^/Mq = 1.0 x 10^ (upper), 2.4 x 10^ 
(heavy, with 95% confidence bands) and 3.8 x 10^ (lower). The velocity dispersion data from 
which these curves were derived extends only to 25"; hence, the confidence bands become very 
wide at large radii. 
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